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STUDIES  OF  ATMOSPHERIC  WAVES  AND  TURBULENCE  PROCESSES 
USING  NUMERICAL  MODELING  AND  ATMOSPHERIC  MEASUREMENTS 

I.  Project  Objectives 

Our  objectives  for  the  original  AFRL  project  were  many.  A  complete  description  of  those  eoals  is 
provided  m  the  last  R&D  status  report  covering  the  period  September  1998  through  Lptember  1999  The 
najowed  scientific  goals  of  the  descoped  research  effort,  extending  from  March  1998  to  December'  1999 

1.  Perform  numerical  simulations  and  analysis  of  Kelvin-Helmholtz  (KH)  shear  instability  assess 
turbulence  structure  and  its  impacts  on  stable  stratification  and  shear;  and 

simulations  to  address  the  2nd  major  source  of  shear  turbulence  in 

II.  Progress 

Significant  advance  were  made  in  both  research  foci  under  this  contract.  The  first  task  resulted  in 
the  most  complete  simulations  of  KH  instability  and  turbulence  generation  in  existence  ThfflrS 
application  of  these  results  was  in  addressing  the  anisotropic  nature  of  stratospheric  turbulence  The 
egree  of  anisotropy  in  stratified  flow  has  been  evaluated  as  a  function  of  depth  through  the  shear  layer 
^d  ^  a  function  of  the  Reynolds  number  (Re)  for  two  values  of  the  RichardL  number  (Ri)  We  SS 
stratification  dominates  the  dynamics  in  the  edge  regions  of  shear  layers.  This  tendency  is  grrater  at 
higher  values  of  Ri  and  persists  for  all  values  of  Re.  The  dynamics  in  these  edge  regions  are  chafacterized 
by  entrainment  with  ambient  fluid  outside  the  turbulent  layer.  As  such  the  value  of  the  effectwI  S 
e^olds  number  is  lower  than  m  the  turbulent  interior  of  the  layer,  and  the  dynamics  are  dominated  by 

coherent  vortex  tutes  onented  pnmarily  with  the  streamwise  direction  near  the  layer’s  edge  and  oriented 
With  the  spanwise  direction  at  the  layer’s  edge.  ^  ®  ana  oriented 

Spectral  slopes  consistent  with  predictions  by  Bolgiano  (1959)  for  stratified  turbulence  are 

of  the  shear-flow  evolution.  However,  at  later  times  when  turbulence 
has  significanUy  eroded  the  primaiy  instability  mode  (i.e.,  at  a  time  when  one  might  expect  better 
agreement  with  Kolmogorov’s  1941  predictions  than  Bolgiano’s  1959  predictions),  the  fpectrchartler 
of  stratosphenc  turbulence  exhibits  a  slope  between  that  predicted  by  Bolgiano  (-1 1/5)  and^Kolmogorov 

cn  I  /’f  velocity  field,  but  which  remains  consistent  with  Bolgiano's  -7/5  prediction  for  the 

spectral  slope  of  the  temperature  field.  Efforts  near  the  end  of  the  contract  Lgan  to  address  the  structure 

functions  accompanying  sheared  turbulence  and  provided  additional  evidence  of  departures  from 

expectations  of  isotropic  homogeneous  turbulence.  Direct  measurements  of  the  turbulence  inner  scale^and 

comparison  with  its  theoretical  dependence  on  energy  dissipation  rate  nevertheless  demonstrated  that  the 
theoretical  relationship  is  highly  accurate.  ciuGusiraiea  mat  me 

Initial  studies  of  wave  breaking  under  the  second  research  task  were  also  completed  but  have  not 
been  evaluated  m  the  same  detail.  Because  of  limited  resources,  our  focus  shifted  from  an  anelastic  code 
development  to  Boussinesq  simulations  for  gravity  waves  having  amplitudes  both  above  and  below  a 
convective  instability  threshold.  Simulations  at  amplitudes  below  L  cLvective  threshold  s"pportedThe 
h„ear  s,ab.h.y  analysis  of  Iximbard  and  Riley  (1995),  suggesting  that  such  motions  can  leaVtZs«lt 

antic«letoarrC^ 
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Some  time  was  also  spent  conducting  preliminary  2D  tests  so  that  more  expensive  3D  simulations 
would  employ  computer  and  personal  resources  optimally.  Some  preliminary  3D  cases  were  run  at  lower 
resolution  to  understand  the  morphology  of  the  gravity  wave  as  it  began  to  break  down  and  generate 
turbulence.  Once  an  understanding  of  the  resolution  needs  was  in  hand,  our  first  "high-Re"  solution  was 
computed  and  has  recently  been  completed.  As  with  earlier  work,  we  observe  the  development  of 
intertwined  hairpin  vortices  at  the  edges  of  the  unstable  regions.  Propagation  of  the  wave  during 
turbulence  evolution  results  in  a  trail  of  mixed  fluid;  unlike  the  turbulent  shear  layer  which  concentrates 
mixing  in  a  fixed  layer. 

III.  Invention  Disclosure  and  Patents  Granted 

None. 

IV.  Property/Equipment 

Three  hard  drives  and  a  memory  kit  for  the  origin  and  onyx  systems  were  purchased.  These  items 
were  placed  on  the  government  inventory  schedule  12/10/01.  The  plant  clearance  case  number  is 
S0602A-1311. 

V.  Personnel  Changes,  Administrative  Actions,  Conferences 

The  supplement  allowed  for  maintenance  of  a  skeleton  presence  to  continue  the  simulations  so  the 
valuable  DoD  supercomputer  resources  we  have  obtained  (450,000  T3E  this  year  alone)  were  not  wasted. 
Dr.  Joe  Weme  presented  work  at  the  European  Geophysical  Society  and  the  Monterey  DoD  HPCMO 
meeting. 

VI.  Travel 

Travel  to  Den  Haag,  Netherlands  and  Monterey,  CA  were  required  for  the  conferences  attended  by 
Dr.  Joe  Werae. 
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